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ABSTRACT 



There is considerable discrepancy between the amount of X-ray absorption and that inferred from optical (rest frame UV) as measured 
along gamma-ray burst (GRB) sight lines, with the former being typically an order of magnitude higher than what would be expected 
from the measurement of neutral element species via optical absorption line spectroscopy. We explore this "missing gas problem" 
by using X-ray and optical measurements in a sample of 29 z = 0.7 — 6.3 GRBs from both spectroscopic data and the afterglow 
broadband spectral energy distributions. The low ionisation species detected in the UV are associated with the neutral interstellar 
medium in the GRB host galaxy, while soft X-ray absorption, which is weakly dependent on the ionisation state of the gas, provides 
a probe of the total column of gas along the sight line. After careful consideration of any systematic effects, we find that the neutral 
gas consists of <J 10% of the total gas, and this limit decreases with the more ionised that the X-ray absorbing gas is, which in our 
spectral fits is assumed to be neutral. Only a very small fraction of this ionised gas, however, is detected in UV absorption lines 
with ionisation potentials up to ~ 200 eV (i.e. Si IV, C IV, N V, O VI), which leaves us to postulate that the X-ray excess is due to 
ultra-highly-ionised, dense gas in the GRB vicinity. 
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1. Introduction 

The vast luminosities emitted by quasi stellar objects (QSOs) 
have long been recognised as providing a powerful tool with 
which to study the environmental properties of high-redshift 
galaxies, and ever since th e discovery of gamm a-ray burst (GRB) 
afterglows (GRB 970228: Ivan Paradiis et aflll997l) GRBs have 
also been increasingly used to probe distant, otherwise inac- 
cessi ble environments through their highly luminous emission 
(e.g. iKawai et all l2006t jgereer et all 12006b: iFvnbo et all l2006t 
Ruiz-Vel asco et all 120071: iGreiner et all l2009t) . Recent years 
have seen an increase in the abundance of early-time, broadband 
photometric and spectroscopic data of GRBs, largely owing to 
the success of the rapid-respon se, multi-wavelength GRB mis- 
sion, Swift dGehrels et alll2004h . and this has revolutionised our 
capabilities to study in detail the GRB afterglow and the proper- 
ties of its host galaxy. GRB observations are n ow reaching fur- 
ther distances than QSO s (e.g. GRB 090423: ISaivaterra et all 
2009; iTanvir et"aQ 120091) . and their transient nature allows the 
host galaxy to be studied at later epochs free of the blinding GRB 
afterglow (e.g. iFruchter et ail 120061 ISavaeho et all \2Q0% . 

There is now a wel l established connection b etween long du- 
ration GRBs (> 2 s; iKouveliotou et all 1 1993b . which are the 
focus of this work, and the formation of massive stars (e.g. 
GRB980425: Galama et al. 1999, GRB 030329: Hjorth et al. 
2003; Stanek et al. 2003). Due to the short life times of mas- 
sive stars, GRB afterglow observations therefore probe the gas 
and dust within areas of star formation, in addition to the inter- 
stellar medium (ISM) of young, distant galaxies, as well as the 
intergalactic matter (IGM) intervening these lines of sight. 



From the increasing sample of GRB spectroscopic observa- 
tions that cover the rest frame Lya feature at 1215 A, it is be- 
coming clear that GRB host galaxies have large column den- 
sities of cold neutral gas (T<; 10 3 K), with the large fraction of 
GRB hosts containing a damped-Lya (DLA) system (log A^hi > 
20.3 cm" 2 ) or sub-DLA (19.0 < log N m < 20.3 cm" 2 ). The 
survival of certain species, such as Mg I, and time varying Fe II 
and Ni II fine-structure lines place this neutra l gas component 
at a f ew hundred parse cs fro m the GRB (e . g. IProchaska et all 
2006; IVreeswiiketall l2007t iLedoux et all l2009t) . within the 
ISM of the host galaxy. 

Absorption by molecular hydrogen in the UV rest frame may 
provide a probe to the star-forming reg ion closer in to th e GRB, 
where H2 would be expected to exist (Tumlinson e t all 1200' 
Nevertheless, de s pite ex erted efforts (e.g. lVreeswijk et afl . 1200' 
Tumlinson et afj, 120071) . molecular hydrogen absorption has, 
as of yet, o nly been clearly de tected in the spectrum of 
GRB 080607 (IProchaska et allboogh 1 . and even then the molec- 
ular hydrogen may have originated from another molecular 
cloud along the line-of-sight distinct from the star-forming re- 
gion surrounding the GRB. However, the number of GRBs 
with adequate spectral resolution and coverage to detect the H2 
Lyman-Werner absorption lines continues to be small, and H2 
absorption lines ste mming from the G RB star-forming region 
may still be detected ( Tu mlinson et al l l2007h . 

For high-z GRBs (z ;> 2.0), the circumburst environment 
may be further probed by the detection of absorption lines pro- 
duced by highly-ionised species (i.e. ionisation potentials be- 



34; 



1 Fvnbo et al. (2006) reported a tentative detection of H 2 in the spec- 
trum of GRB 060206. However, iTumlinson et all j2007l) later argued 
the claimed H2 absorption lines to be Lya forest lines. 
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tween 50-150 eV) present in the hot gas (T~ 10 4 K) surround- 
ing the GRB. Although these highly ionised absorption lines 
may also originate from hot gas in the halo of the galaxy, mea- 
surements of O VI, C IV, Si IV and N V in a handful of GRB 
afterglow spectra hav e been used to probe the gas close in to t he 
GRB progenitor star dFox et all 12001 iProchaska et alll2008l) 

In contrast to the specific regions of gas that can be iden- 
tified from UV spectra, X-ray spectroscopic observations pro- 
vide measurements of the total column density of gas along the 
line-of-sight, probing both the cold neutral gas as well as the 
warm ionised regions. Soft X-rays with energies < 0.8 keV 
are absorbed by medium-weight metals along the line-of-sight, 
which is predominantly in the form of oxygen and to a lesser ex- 
tent carbon and nitrogen. Furthermore, the cross-section of oxy- 
gen remains relatively unc hanged regardless of ionisation state 
dVerner & Yakovlevi 1 1995b . making it a good proxy for the host 
galaxy total oxygen column density. 

In addition to the conditions of the host galaxy gas along 
the line-of-sight, the host dust abundance and extinction prop- 
erties can be explored from dust depletion patterns inf erred 
from relative metal abundances (e.g.lSavaglio & Fall, 2004), and 
also from the analysis of broadband infrar ed (IR) through to 
UV spectral energy distributions (SEP) (e . g. | Galama & W ii ers , 
200lUStratta et al.Ll2004tlKann et alll200alstarling et al.Ll2007t 
Schad v et all 120071 l2010al) . Soft X -rays are also extinguished 
and scattered by dust, although by much smaller amounts such 
that absorption by mediu m weight metals is the dominant sourc e 
of attenuation (e.g. Weingart ner & DraineLl200ll:lDrainell2003h . 

Due to the complex structure of the absorbing gas and the 
uncertainties in the location of the various attenuating gas and 
dust components, a combined analysis of the UV and the soft 
X-ray energy bands gives a more comprehensive understanding 
of the properties of the circum- and interstellar medium of GRB 
host galaxies. For example, the order of magnitude larger equiv- 
alent neutral hydrogen column densities typically derived from 
soft X-ray absorption 2 , A^h,x, in comparison to the host galaxy 
neutral hydrogen colum n densities measur ed from damped Ly- 
a absorption, TVhi, led lWatson et al.l d2007l) to conclude that the 
former probed a significant column of ionised gas that is trans- 
parent to UV photons. 

In this paper we combine the results from the analysis of op- 
tical and X-ray spectra, and also of broadband afterglow SEDs, 
with the aim of providing greater insight on the ionisation state 
and dust properties of the host galaxy medium along the line- 
of-sight to GRBs. Instead of neutral hydrogen, we use absorp- 
tion lines of singly-ionised metals to trace the neutral gas within 
GRB host galaxies that can then be compared directly with the 
soft X-ray absorption, which measures primarily the oxygen col- 
umn density (neutral and ionised), A^o.x- This removes the need 
for hydrogen column density and host galaxy metallicity mea- 
surements, which are frequently uncertain or unavailable, thus 
increasing our sample size. The results presented here contribute 
towards not only a better understanding on the nature of galaxies 
hosting GRBs, but also on the environmental conditions within 
high-redshift galaxies as a whole. 

In section [2] we provide our sample selection criteria, and 
we describe our data reduction and analysis in section [3] We 
present our results and discuss their implications in section 21 
and our conclusions are summarised in section [5] All quoted er- 
rors throughout the paper ar e I -a. Throughout thi s paper we use 
the cosmic abundances from Aspl und et al.l (120091) . 



Assuming solar abundances from And ers & Grevessel d!989h . 



2. Sample 

The sample selection is based on the requirement that optical 
spectroscopic data are available for the GRB afterglow with ab- 
sorption measurements in at least Zn II, Si II, S II or Fe II with 
which to probe the host galaxy neutral gas. All these low ionisa- 
tion species have ionisation potentials (IPs) slightly larger than 
hydrogen, and all are frequently measured in GRB optical spec- 
tra. Our initial sample amounted 7 pie-Swift and 22 Swift GRBs, 
of which 26 had X-ray afterglow spectroscopic observations (all 
Swift GRBs plus four pie-Swift GRBs), and 23 had optical/X- 
ray SED-fit host galaxy visual extinction measurements, Ay, (all 
pie-Swift and 16 Swift GRBs) available either from the literature, 
or with sufficient UV through to infrared (IR) data for an SED 
analysis of our own to be possible. Due to the degeneracy that 
exists between spectral index and extinction, we only consider 
those Ay values derived from broadband optical to X-ray SED 
modelling. As a consequence of our selection criteria, whereby 
at least one metal line is detected in the optical spectra, a spec- 
troscopic redshift measurement is available for all GRBs in our 
sample. 

Our requirement for an optical afterglow spectrum intro- 
duces a selection effect against dim or highly extinguished 
GRBs, either from dust in the host galaxy, Galactic dust, or 
dust within intervening systems. GRB host galaxy ISM prop- 
erties are independent of the Galactic or foreground dust ex- 
tinction, and selection effects biased against GRBs with large 
foreground extinction should, therefore, not affect our overall 
results. Similarly, selection effects against redshift ranges that 
have few prominent absorption lines in the observer frame op- 
tical bandpass (the 'redshift desert'; z ~ 1 — 1.5) are likely in- 
dependent of the GRB host environmental conditions, and again 
should not affect our overall findings. 

The dominant selection effect that we need worry about is 
the bias against GRBs highly extinguished within their host 
galaxy. Recent studies suggest that between 20-30% of GRBs 
are undetected in the optical wav e length range as a result 
of dust extinction (e.g. ICenko et all 120091: iFvnbo et ail [2009; 
iPerlev et all 120091: iGreiner et all 1201 0). and there is some in- 
dication that these more extinguished GRBs reside in more 
massive and evolved galaxies than the hosts to optically bright 
GRBs. It should, therefore, be kept in mind that the results pre- 
sented in this paper apply to GRBs with low host galaxy extinc- 
tion along the line-of-sight (Ay <^ 1.0). In the case of selection 
effects against dim GRBs, our results will be affected if low lu- 
minosity GRBs belong to a different cla ss of GRB to their high 
lumi nosity counterparts , as su ggested by Liang & Zhang ( 2006) 
and Nardini et al . (2006t l2008l) (although see lOates et alll2009l) . 
and this possibility should, therefore, be kept in mind. We, there- 
fore, stress that the results presented in this paper may not apply 
to low luminosity and/or highly extinguished GRBs. 

3. Data analysis and reduction 

All optical spectroscopic column density measurements used in 
this paper are listed in Table [TJ and have mostly (~ 75%) been 
taken from the literature, in which case the corresponding ref- 
erence is also given. In a few cases we re-computed the col- 
umn densities ourselves using the published equivalent widths 
(EW) and applying standard techniques. For the 26 GRBs with 
X-ray spectroscopic data, we preferentially used the X-ray col- 
umn densities measured from optical to X-ray SED analysis for 
consistency with the Ay vaue that we used. Where this was not 
possible due to the lack of optical/NIR data (five GRBs), we used 
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Table 1. GRB host galaxy column densities for 29 GRBs. Columns list redshift, Galactic column density and visual extinction in the line-of-sight to the GRB, the H I, Si II, S II, 
Fe II and Zn II column density with corresponding references, the oxygen column density measured from the X-ray spectra, No t x, and the best-fit visual extinction, Ay, with 
both with corresponding references given. For GRB 030226 and GRB 090323, two absorbing components were detected within the GBR host galaxy, and the optical column 
densities for both are given where applicable. 
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1 ISavaglio et alj d2003|); 2 IStarling et al.1 d200'ft; 3 iHjorth et al.1 d2003l); 4 iMasetti et al.1 d2003l); 5 This w ork; 6 ISavaglio & Fail d2004l) ; 7 IShin et al.l d2006l); 8 IVreeswiik et alJd2004b: 9 IWatsonetal 
2006): 1U lEvans et al.1 d2009h: 11 iBerger et alld2006f): " iLedoux et all d2009h: 13 ISchadv etafl d2010bl); 14 [Prochaska et al. [ 2007ah; 15 lTotanil d2006h:" Tg iKawai et al.1 d2006h: 1 ^TPiranomonte et al 



2008) 



18 )Fvnbo etalj d2006l); 19 [Th6ne et al.1 (2008a) 
t alJd2010h : 27 iPerlev et all d2008a 



2009); ^ D lGreineret; 



20 lAokietalJ d20 09): z 
; 2& lD'Elia et alj d2009l) ; 



Vreeswiik et al. ( 2007); " Chary et a ll d2007l); iPrice et al.1 d2007p: ^ |Th6ne et alj d2008b[) ; a Elfasdottir et al 

309); '" 



2a lD'Eliaeta l. (2009); 30 iProchaska et all d200' 



Starling ct 



:007); nil 
alJ d2009l)~ 



' iMcBreen etaiTd2010l) ; •" iRau et al 



20ld) 

* Where two host galaxy absorption systems are reported, we used the combined column densities of the two systems in our analysis. 
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X-ray absorptio n values measured f rom the X-ray data alone, 
which we took from Evans et al ] (120091) . 

Both in our own X-ray analysis and for the res ults taken from 
the literature, the X-ray fitting package XSPEC (lArnaudl fl996) 
was used to model the X-ray spectra with an absorbed power- 
law (phabs*pow), where solar abundances were assumed. This 
provided an equivalent neutral hydrogen column density, A^h,x> 
which we converted back into an oxygen column density, Nq ; x> 
and also list in Table Q] with corresponding references where rel- 
evant 3 . We note that in using the photoelectric absorption model 
phabs we assume that the soft X-ray absorbing gas is neutral , 
which may not be an accurate approximation (e.g IWatson et al.L 
120071: iProchaska et all [2008 ), in which case the column of ab- 
sorbing material would be underestimated. The best-fit A^h,x, 
and consequently also A^o,x listed in Table Q] are thus effective 
lower limits. Furthermore, for increasing redshift, the absorption 
from oxygen is increasingly shifted out of the bandpass such that 
at z > 2.5 the majority of the soft X-ray absorption is caused by 
other heavier elements. In this case our oxygen column density 
measurements are indirectly determined assuming solar relative 
metal abundances, which is an assumption typically made in X- 
ray spectroscopic analysis. However, the implications of such an 
assumption apply as much to A^h.x as they do to iVb.x, but in the 
latter case a lot of the uncertainty corresponding to the metallic- 
ity of the environment for the large fraction of GRBs at z > 2.5 
is removed. 

In the final column of Table Q] we give the best-fit host 
galaxy visual extinction along the line-of-sight to the GRB, Ay, 
with relevant references provided when taken from the literature. 
Finally, the GRB redshift, and the G alactic soft X- r ay ab sorp- 
tion and visua l extin ction taken from Kalb erla et alJ (120051) and 
Schle geTet alj (l998), respectively, is also given in Table [TJ 



3.1. SED Analysis 

For consistency within our sample we require that the host 
galaxy visual extinction and soft X-ra y absorption are mea - 
sured in the same w ay fo r each GRB. In Schad y et alj d2010al) . 
iGreiner et alj d2010h and Starling et al] d2007l) . host extinction 
and absorption values were measured from the same model fits 
to the GRB SEDs, which we describe in the next section. We, 
therefore, consider the results from these papers comparable to 
each other. For six GRBs in our sample not analysed in any of 
the papers listed above, sufficient optical/NIR and X-ray data 
was available to produce a broadband SED. In these cases we 
therefore pe rformed our own SED analysis following the method 
described in lSchadv et alj d2010al) . 

GRB 020405 and GRB 030226 had available Chandra 
data, which we reduced using the standard procedure. The 
optical data for these SEDs was taken from the literature 
(GRB 020405; Masetti et al. 2003, GRB 030226; Klose et al. 
2004). The remaining four GRBs (GRB 050401 GRB 090323 
GRB 90328, GRB090 926A) had S wift X RT .Burrows et all 
l2005ah and UVOT dRoming et all [2005) data, which were 
taken from the UK Swift data archive 4 , as well as g'r'i'z' and 
JHK da ta taken with the GR OND, a simultaneous 7-channel 
imager dGreiner et all |2008) mounted at the 2.2m MPI/ESO 
telescope at La Silla (Chile). All XRT data were in photon 
counting (PC) mode and were reduced with the XRTPIPELINE 



tool (vO.12.4 ). Spectra were ex tracted following the method 



described in Evan s et al ] (I2009h . avoiding time intervals that 



3 By default XSPEC uses the abundan ces from|An ders & Grevesse 
( 1989), which are consistent with those of Aspl und et al.1 d20(M) within 
errors. Variations in Nu,x and iVo,x measurements resulting from the 
use of different abundances are, therefore, negligible. 



showed evidence of spectral evolution, such as du ring the early- 
time steep decay phas e of the X-ray light c urve dNousek et all 
2006), during flares dFalcone et all 120071) . Both UVOT and 
GROND dat a were reduced fo llowi ng the standard procedures 
described in iPoole et ail d2008l) and iKupcii- Yoldas et alJ (2008) 
respectively. 

The SEDs were fitted within XSPEC (vl2.6.0) with a Galactic 
and host absorbed and extinguished power law, allowing for a 
cooling break between the X-ray and optical e nergy bands (with 
change in spectral slope fixed to A/3 = 0.5; ISari et all fl998) 
when it provided a significantly improved fit. Neutral hydrogen 
absorption from the IGM was als o taken into a ccount in our 
model using the prescription from M adaul d 19951) . Solar abun- 
dances are assumed in the soft X-ray absorption measurements, 
and the host visual extinction was measured from the best-fit to 
the data between the Small or Large Magellanic Cloud (SMC 
and LMC, respectively), or the Milky Way mean extinction laws 
dPelll992h . 

3.2. Optical Spectroscopy Analysis 

For the ~ 25% of the sample where we performed our 
own optical spectroscopic analysis, column densities of met- 
als wer e deriv ed by the curve of growth (COG) analysis 
(Spitzer, 1978 ), or if appropriate , the apparent optical depth 
(AOD)( Savage & Sembachl 1199 lh method. The COG provides 
reliable results through the linear approximation when the equiv- 
alent width of a line is EW < 0. 1 A and the effective Doppler pa- 
rameter is b > 20 km s _1 . For stronger lines, the COG requires 
the detection of at least several lines with different oscillator 
strengths of the same ion. We, therefore, used the AOD method 
for moderate to low EWs, provided that the effective Doppler 
param eter is not very low (6 > 10 km s _1 ; IProchaska et all 
l2007ah . 



4. Results and Discussion 

4. 1. Dust Within the GRB Host Galaxy 

Iron is a highly refractory element and its relative abundance 
with non-refractory, iron-peak elements acts as a good tracer 
of the abundance of dust particles present within the region of 
absorbing gas. The column density of singly ionised iron and 
zinc are often used for this purpose since both are frequently 
measured in UV spectra and they both trace neutral gas. It has 
been suggested that a relative zinc-to-iron abundance of up to 
[Fe/Zn] 5 w -0.3 dex could arise in low-metallicity environments 
due to differences in the n ucleosynthesis processes from solar 
(IProchaska & WolfeL l2002h . However, any zinc-to-iron relative 
abundances lower than this are likely to be indicative of dust de- 
peletion. 

In Fig. Q] we plot the iron-to-zinc relative abundance [Fe/Zn] 
as a function of the Zn column density, Nz n u, for 18 GRBs in 
our sample with Fe II and Zn II column density measurements, 
where at least one is well constrained (i.e. not an upper or lower 
limit). Of these 18 GRBs, five have only an upper limit on the 
Zn II column density and these are plotted with smaller sym- 
bols, and also with l eft- and upward pointing arrows. A sample 
of QSOs taken from Savaglio] d2006l) are also plotted with open 



http://www.swift.ac.uk/swift_portal/ 
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Fig. 1. Iron-to-zinc relative abundance, [Fe/Zn], as a function of 
the logarithmic Zn II column density. GRB host galaxy measure- 
ment s are plot t ed as solid circles and a sample of QSOs taken 
from Savaglio (2006) (open circles) are plotted for comparison. 
The dashed line shows the ordinary least-squares (OLS) bisec- 
tor fit to the combined GRB and QSO data, which has a power 
law index of -0.7. GRBs with limits on either Fe II of Zn II are 
plotted in grey and with smaller symbols. 



circles for comparison, and the dashed line is the line of best-fit 
to the combined GRB and QSO sample. Given the increasing 
evidence indicating that the bulk of neutral gas within GRB host 
galaxies probed by optical afterglow spectra lies at a few hun- 
dred parsecs from the GRB, the similarity in trend of smaller 
iron-to-zinc relative abundances for larger Zn II column densi- 
ties between QSOs and GRBs is not so surprising. A Spearman 
rank test for the 13 GRBs with well measured Fe II and Zn II 
column density measurements gives a Spearman coefficient of 
p = —0.51 at 90% confidence, and the Spearman rank coeffi- 
cient remains the same when we include the QSO data, although 
the confidence in the coefficient increases, indicating a signifi- 
cant anti-correlation between [Fe/Zn] and N% n u, and thus sug- 
gesting that dust depletion tends to be higher for larger Zn II 
column densities. 

For all 13 GRBs that have well-constrained iron and zinc col- 
umn densities, seven also have measurements of the host galaxy 
visual extinction, Ay, at the 3cr level (see Table [TJ. To test 
whether the correlation between dust depletion and Zn II col- 
umn density indicated by Fig. [T]is in agreement with the best-fit 
Ay values from our SED analysis, we split this subsample of 
seven GRBs into two, depending on their Zn II column density. 
We use a threshold of Nznii = 13.5, which is the mean Zn II col- 
umn density for the GRB subsample, and find that for the set 
with log iVznii < 13.5 and log Nz, n ii > 13.5, the mean best- 
fit Ay values are (Ay) = 0.17 and (Ay) = 0.62, respectively, 
with respective standard deviations of 0.07 and 0.27. Although 
the scatter in Ay is large, it is reassuring that the trend of larger 
dust-depletion for larger Zn II column densities observed in our 
optical spectroscopic analysis is similarly present in our SED 
analysis results. 

This general agreement is better illustrated in Fig. [2] where 
we have plotted the relative iron-to-zinc abundance, [Fe/Zn], 
against visual extinction, Ay, for 16 of the GRBs included in 
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Fig. 2. Iron-to-zinc relative abundance, [Fe/Zn], against visual 
extinction, Ay, for a sample of 16 GRBs. GRBs with limits on 
either [Fe/Zn] or Ay are plotted in grey and with smaller sym- 
bols. 



Fig. Q] that have Ay measurements available. GRBs with lower 
limits in either [Fe/Zn] or Ay are plotted with smaller sym- 
bols, and GRBs with SEDs best-fit with a spectral break be- 
tween the optical and X-ray energy bands are plotted as open 
circles. There is one clear outlier in this GRB sample with a 
large visual extinction of Ay > 1.0 but small [Fe/Zn]. This data 
point corresponds to GRB 070802, and an abundance of dust 
along its line-of-sight, as suggested by it's large afterglow ex- 
tinction, may thus no longer make our assumption of negligible 
zinc depletion valid. Further evidence of the line-of-sight to this 
GRB having significantly amounts of dust is in the clear detec- 
tion of a 2 175 A dust extinction feature in its optical spectrum 
dKriihler et all l2008t lElfasdottir et all l2009) 6 . which is thought 
to be more prominent in du sty environments (e.g. iNoll et all 
12007 : ProchaskaetaJi |2009). Although along numerous lines- 
of-sight the fraction of zinc in grains is ge nerally negligible, it 
can b e as high as 20% in the cool disc ISM (Savage & Sembach, 
1996), where the dust abundance is higher. The inconsistencies 
between [Fe/Zn] and Ay for GRB 070802 may therefore be due 
to a significant amount of both Fe and Zn dust depletion, causing 
the relative iron-to-zinc metal abundance to only be a little lower 
than solar abundances. 

4.2. Gas lonisation State 
4.2.1. Total vs. Neutral Gas 

To trace the neutral gas within a GRB host galaxy, we prefer- 
entially use the ion Zn II or S II, which if not available is sub- 
stituted by Si II or otherwise Fe II. This is due to the inherent 
uncertainties present in the dust depletion of the refractory ele- 
ments Si II, and especially Fe II. When column densities for both 
Zn II and S II are available, the one with the smaller error bars 
is used. When only Si II or Fe II absorption measurements are 
available, a dust correction is applied to the measured column 

6 GRB 08060 7 als o had a well-detected 2175 A feature 
dProchaska et all 120091) . However, only lower limits on the zinc 
and iron column densities are available for this GRB. 
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Fig. 3. Host galaxy neutral gas column density, iV n tr> against 
host galaxy total gas column density, A?n,o> along the line-of- 
sight to a sample of 26 GRBs. For each GRB, N ntT is derived 
from either the column density of Zn II (red circles), S II (or- 
ange squares), Si II (green stars) or Fe II (blue triangles), where 
a correction for dust depletion is applied to A^.si and A^Fe- 
Smaller data points correspond to those taken from low- or mid- 
resolution spectra, and larger data points are taken from high- 
resolution spectra (R > 10, 000). The dashed line corresponds 
to where N ntr is equal to A^o 



density. The general consistency between the dust depletion in- 
ferred from iron-to-zinc relative abundances and the SED best-fit 
visual extinction shown in the previous subsection gives confi- 
dence to the assumption that zinc dust depletion is typically neg- 
ligible. We would expect to see a discrepancy between [Fe/Zn] 
and Ay where this assumption no longer applies (see Fig. [2]). 
To correct for the dust depletion of iron and silicon, we there- 
fore fitted a linear function to logiVsui and logA^oH against 
log iVz n n, respectively, for the combined data set of GRBs and 
QSOs, and used these lines of best-fit to correct for dust deple- 
tion in our iVsni and A^eii measurements. We then corrected 
for cosmic abundance variances by normalising al l our measure- 
ments to each other using the abundances from lAsplund et al.l 
(2009). We denote the normalised column densities of oxygen, 
zinc, silicon, sulphur and iron by Nn,o> Nn,Zm Nn.Si, Nn,s 
and Njf p e , respectively. 

In Fig. [3] we plot the neutral gas column density, N n % T , rep- 
resented by either A^.zn (red circles), A^n.s (orange squares), 
-An, si (green stars) or A^n.Fc (blue triangles), against N^ t o, 
which represents the total column of gas (neutral and partially 
ionised), where both A^si an d A^n.Fc have been corrected for 
dust depletion. We have also differentiated between the reso- 
lution of the optical spectra used for each GRB observation, 
where those measurements taken from high-resolution spectra 
(R > 10, 000) are plotted with larger symbols, and measure- 
ments from low- or mid-resolution spectra are have smaller sym- 
bols. Just under two thirds of our sample have low- or mid- 
resolution spectral data, and this may have consequences for the 
precision of the column density measurements in these cases, 
which we discuss further below. The dashed line in Fig. [3] cor- 
responds to where the column density of neutral gas is equal to 
the total column density of gas as probed by our X-ray absorp- 



tion measurements. From this figure it is clear that for all GRBs, 
N?s o is larger than the column density of neutral gas, AT ntr , by 
typically around an order of magnitude. 

4.2.2. Systematic Effects 

The neutral gas column density in two thirds of the GRBs plot- 
ted in Fig [3] was measured from Fe II and Zn II absorption, 
and in these cases at least some of the excess absorption mea- 
sured in X-rays could be due to GRB hosts having a-element 
overabundances, a-elements are produced within the sequence 
of He, such as O, Mg, Si and S, and trace massive star forma- 
tion. An over-abundance of a-elements should, therefore be ex- 
pected in young , star-fo rming galaxies such as GRB hosts (e.g. 
ProchaskaeialJ,|2007b). However, such overabundances would 
not account for the order of magnitude differences between X- 
ray and optical absorption measurements, nor would it apply to 
the third of the sample with neutral gas column densities traced 
by Si II absorption. We, therefore, do not consider an abundance 
in a-elements within GRB host galaxies to be the dominant 
cause of the larger values of A^o measurements compared to 
AW 

Alternatively, saturation and lack of resolution in the opti- 
cal spectral data could lead to the metal column densities be- 
ing underestimated (e.g lProchaskal I2006I) . However, this is not 
a problem when high-resolution spectra are available, which is 
the case for ~ 40% of our sample. Optical column density mea- 
surements for this fraction of the sample are plotted with larger 
symbols in Fig. [3] and most of these fall well below the dashed 
line, in the region of space where the total gas column density 
is significantly larger than the neutral gas column density. For 
the remaining ~ 60% of the sample, the error on the neutral gas 
column density measured from low- and mid-resolution optical 
spectra may be as large as 0.6 dex, although this would be a 
highly pessimistic case. However, even in this unlikely situation, 
over half of the mid- and low-resolution data points would still 
have a total gas column density at least an order of magnitude 
larger than the neutral gas column density. We therefore con- 
clude that although the resolution of the optical spectroscopic 
data available may cause the neutral gas column density to be 
under-predicted, for the majority of the data points in Fig. [3] the 
neutral gas column density would remain significantly smaller 
than the metal column density implied by the X-ray data. 

A clear possibility that we need to address is that of un- 
accounted intervening absorbing systems in our line-of-sight, 
which would result in the GRB host soft X-ray absorption be- 
ing over-estimated. In optical spectroscopy the location of any 
intervening absorption systems are clearly determined, and thus 
the absorption from such systems is easily separated from ab- 
sorption within the GRB host galaxy. However, the resolution of 
our X-ray spectra is insufficient for the contribution from several 
absorbing systems on the GRB afterglow soft X-ray absorption 
to be disentangled. It should be said that such strong intervening 
absorbing systems are rare, and the probability of there being 
an absorption system in the line-of-sight to a GRB of sufficient 
optical depth to contribute a mea surable amount of soft X-ray 
absorption is, t herefor e, small ( O' Flaherty & Jakobsen, 1997; 
ICampana et all l2006bl) . In IWatson et al.l ([2007) this issue was 
looked at in greater detail, and they estimated that the absorp- 
tion resulting from any intervening system would need to be two 
orders of magnitude larger than the absorption measured thus 
far in optical spectra; specifically than the absorption measured 
from Mg II sys tems, which are typi cally the largest of absorbers. 
Further to this, Wat son et al.l d2007h noted the soft X-ray absorp- 
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tion in the line-of-sight to blazars to be small relative to that of 
GRBs. This is despite the two types of sources having compara- 
ble numbers of intervening Mg II absorption systems along the 
line-of-sight, thus further supporting an intrinsic absorption in- 
terpretation. 

A greater worry in our X-ray measurements is that of the 
Galactic absorption, which i n the case of our G RB sample can 
be uncertain by up to 15% (Kalb erla et all 120051) . Any under- 
estimate of the Galactic soft X-ray absorption would result in the 
host galaxy soft X-ray absorption being o ver-estimated by a fac- 
tor of (l + z) 2 ' 6 dGalama & W iiers. 2001). In order to investigate 
whether such a factor could account for the large values of -/Vn,o 
observed in Fig. [3] we re-determined the host galaxy X-ray ab- 
sorption for each GRB plotted in the figure, but this time setting 
the Galactic absorption along the line-of-sight to its 3c upper 
limit. We found the re-calculated host galaxy X-ray absorption 
values to have typically decreased by a little less than 50 %, thus 
remaining a factor of five or so larger than iV n t r . Any uncer- 
tainties in our Galactic absorption corrections therefore have no 
significant impact on our primary result that -/Vn.o is a factor of 
several larger than A ntr . 

We therefore conclude that the offset between An,o and 
iV n tr within GRB host galaxies is predominantly due to X-ray 
absorption measurements probing a larger column of gas than 
the UV low ion absorption line measurements. This larger col- 
umn of gas must be ionised, which is not traced by the low ion- 
isation species that we have analysed, but which does absorb 
photons in the soft X-ray energy bands. Dust extinction of soft 
X-rays may also contribute to the large X-ray inferred column 
densities. However, due to the effect of self-shielding, the frac- 
tion of soft X-ray photons lost to extinction will be significantly 
below the frac tion absorbed by oxygen, carbon and nitrogen (e.g. 
iDrainel l200l . Furthermore, given that the soft X-ray absorp- 
tion column densities were measured from a neutral absorption 
model fit to the X-ray spectra, a large ionised gas absorption 
component would turn the X-ray column density measurements 
into lower limits, as already indicated in section [3] The order of 
magnitude difference between the column density of neutral gas 
and the total gas column density measured from the X-ray data 
would thus suggest that, at least along the line-of-sight to the 
GRB, ~ 90 % of the gas within GRB host galaxies is ionised. 

In at least some GRBs most of th e neutral gas lies at a 
few hundred parsecs fr om the GRB (e.g. Prochaska et al., 2006; 
IVreeswiik et alll2007l) . suggesting that little neutral gas remains 
within the star forming region in which the GRB is likely to be 
embedded. It may therefore not be so surprising that the ratio 
between ionised-to-neutral gas is so high. As an illustrative ex- 
ample, if we assume a GRB to lie within a dense star-forming 
region of radius 20 pc then a hydrogen density p w 50 cm -3 
within this region would be needed to account for our soft X-ray 
observations, which is not unreasonable. Such a scenario pro- 
vides a simple explanation to the larger column densities of soft 
X-ray absorbing gas compared to the neutral gas column den- 
sities. It is also consistent with the association of long GRBs 
with massive star formation, whereby GRBs are expected to lie 
within dense star-forming regions of the host galaxy, as well as 
being compatible with the neutral gas lying several hundred par- 
secs from the GRB, as indicated by optical spectroscopy of some 
GRB after glows. 

In both IWatson et all (120071) and iProchaska etal] (120081) the 
quoted distances from the GRB out to which gas could be 
ionised were smaller than the radius of 20 pc that we consider 
above (3 pc and 10 pc, respectively). However, these smaller 
distances were based on the assumption that the GRB was the 
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Fig. 4. Logarithmic host galaxy column density of the highly 
ionised atoms C IV, Si IV, N V and O VI against the total log- 
arithmic column densities of C, Si, N and O in the top, second, 
third and bottom panels, respectively. The total column densities 
are derived from X-ray observations, denoted by Nx- The C IV, 
Si IV , N V and O VI m easurements plotted as circles are all taken 
from Fox et all (I2008I) the N v data plotted as squares are from 
IProchaska et all d2008l) . and data taken from lD'Elia et alj d2010h 
are plotted as triangles. In all four panels, the dashed line cor- 
responds to where the normalised soft X-ray column density is 
equal to the column density of the corresponding highly ionised 
atom. 



sole source of ionising photons. When considering the contri- 
bution from other massive stars in the neighbourhood, the ra- 
dius of th e ionisation fron t could increase out to > 30 pc from 
the GRB (IProchaska et all I2008L although IWatson et all (I2007I) 
argued there only to be a small increase in size of the ionised 
region). Nevertheless, even in the case where the size of the 
ionised region is within 3-10 pc of the GRB, a particle density 
of p w 100 — 250 cm -3 , which remains within reason, would 
similarly satisfy our observations. 



4.2.3. Highly Ionised Gas 

To investigate this scenario further, we take a look at the col- 
umn densities of metals in a higher ionisation state than has 
been looked at so far. Whereas low ion elements will only ex- 
ist far from the GRB, such as Zn II, Si II and Fe II, more highly 
ionised metals, such as C IV, Si IV, N V and O VI could sur- 
vive within the circumburst en vironment of the GRB dFox et all 
2008; Prochaska etal., 2008), and may thus make up a large 
fraction of the ionise d gas p robe d by our soft X- r ay ob serva- 
tions. Both iFox et all d2008l) and IProchaska etal] (|2008) used 
the measurements of highly ionised metals to probe the envi- 
ronmental conditions close to the GRB for a combined sample 
of nine GRBs, where the former used the column densities of 
C IV, Si IV, N V and O VI to trace the highly ionised gas, and 
Prochask a et all ((2008) focused on measurements of N V, which 
are typically less saturated. 

To quantify the fraction of gas probed by the soft X- 
ray energy bands that could be in a hi ghly ionis e d stat e se- 
lected the sample of eight GRBs from IFox et all d2008l) and 
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Prochaska et al. (2008) that had soft X-ray absorption measure- 
ments, as well as GRB 090926A , which had high- ion column 
densities measurements reported [D'Elia et al] d2010l) . This sam- 
ple allows us to compare the soft X-ray column densities to the 
highly ionised column densities of N, C, Si and O. To do this 
we normalised the soft X-ray absorption measurements to the 
cosmic abundance of the element in question being compared 
(i.e. either Si, C, N or O). Our results are shown in Fig. [4] 
which compares the logarithmic column density of the highly 
ionised atoms Nqiv, Nsiv, N^v an d Novi against the soft 
X-ray absorption column densities normalised to C, Si, N and O 
cosmic abundances in the top, second, third and bottom panel, 
respectively. The normalised soft X-ray absorption column den- 
sity is denoted by Nx to indicate the column density of C, Si, 
N or O as determined from our soft X-ray absorption measure- 
ments. The colu mn densities of highly ionised atoms taken from 
iFox et all (2008) are p lotted as circles, the N V measurements 
from Prochas ka et al.l (12008) are shown as squares, and those 
data from lD'Elia et al.l (12010) are plotted as triangles. In all four 
panels the dashed line represents where the normalised soft X- 
ray column density measurements are equal to the column den- 
sity of the highly ionised atom represented in each correspond- 
ing panel. 

It is clear from this figure that all data points lie several or- 
ders of magnitude below the dashed lines, indicating that the 
highly ionised gas within GRB hosts only makes up a small frac- 
tion (< 0.01%) of the soft X-ray absorption. Fig. [3] on the other 
hand, indicates that ~ 90% of the gas probed by the soft X-rays 
is ionised. If the large fraction of this gas is in a lower ionised 
state, then there should be a signature of this in the detection of 
strong absorption lines from intermediate ionisation lines with 
IPs between ~ 20 eV and ~ 50 eV (i.e. Si II and Si IV). This 
is, however, not the case. For example, Al III, which has an IP 
of 28 eV, is frequently observed to be weaker than Al II, which 
has an IP of 19 eV. This would, therefore, suggest that most of 
the ionised gas probed by the soft X-rays is in a ultra-ionised 
state, with IPs larger than ~ 200 eV. The signature left by such 
a ultra-ionised gas would lie predominantly in the soft X-ray 
energy range, which is already heavily absorbed by oxygen and 
carbon (both neutral and ionised). Resolving the absorption from 
an ultra-ionised gas is, therefore, beyond the spectral capabilities 
of current fast-response X-ray telescopes a nd will require a fu- 
ture era of X-ray telescopes such as Xenia dPiro et alll2010l) . 

Another possibility for detecting this ultra-ionised gas would 
be in the form of variability in emission lines as the gas cools 
on the order of years (e.g. iPerna et all 120001) . To date, analy- 
sis investigatin g such cooling has only be en done in the case of 
GRB 990712 (Kup cii Yolda? et all 120061) where five epochs of 
observations of the host galaxy taken up to 6 years after the burst 
were used to search for variations in the emission line fluxes. No 
evidence of variability was found, from which the authors were 
able to put a limit on the density of the emitting (or ionised) re- 
gion of n 6 x 10 3 cm~ 3 . This remains consistent with our 
analysis. 



4.2.4. Gas Ionisation Profile 

A number of theoretical studies investigating the effect of a 
GRB on its surrounding environment show the intense radia- 
tion field produced by a GRB to be capable of destroying dust 
out to ~ 10 pc and ionising gas out to te ns of parsecs (e.g. 
IFruchter et alll200ll : lPer na & Lazzati, 2002). The large fraction 
of ionised gas along the line-of-sight to GRBs derived from our 
optical and X-ray absorption analysis may therefore not be that 
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Fig. 5. Ionisation profile of gas within host galaxy of six GRBs 
as traced by nitrogen. Ionisation state of gas is divided into three 
broad components; neutral (N ntr ), highly ionised (Ny), and 
ultra-ionised (-/V u i tr ion)- Highly ionised nitrogen, N V, is mea- 
sured directly from optical spectra, whereas the column density 
of neutral nitrogen is determined from either S II (black circles) 
or Zn II (grey circles), and the ultra-ionised column density of ni- 
trogen gas is determined from X-ray measurements. The down- 
ward pointing triangle represents an upper-limit. 



surprising. What may be of more surprise is the ionisation pro- 
file of the gas, which our analysis suggests has a well detected 
neutral component, a significantly smaller highly ionised com- 
ponent with IPs reaching up to ~ 140 eV, and then possibly a 
very large ultra-ionised component at IPs of > 200 eV. 

As an illustrative example we used the element nitrogen to 
visualise better the ionisation profile of the gas within GRB 
host galaxies along the line-of-sight, which is shown in Fig. [5] 
Nitrogen is detected in a highly ionised form in six GRBs in 
our sample, and the neutral {N ntr ) and ultra-ionised (A^ u i tr ion) 
fraction of nitrogen are derived from the column densities of 
neutral gas as traced by S II (black circles) or Zn II (grey cir- 
cles), and from the gas column density inferred by soft X-ray 
absorption measurements, respectively. These column densities 
are then normalised to the cosmic abundances of nitrogren. An 
upper-limit on the ultra-ionised gas column density, iV u it r ion> m 
Fig (top curve) is represented by a downwards pointing trian- 
gle. 

The ionisation profile shown in Fig. [5] may arise from a 
distribution of gas with a density profile approximating a step- 
function, such that close to the GRB the density of the gas is sig- 
nificantly larger than far from the explo sion sight, as suggested 
in section l4.2.2l|Perna & Lazzatil (f2002) showed that a GRB can 
easily photoionise a dense column of gas (> 10 3 cm -3 ) out to 
a few parsecs. However, at lower densities further away from 
the GRB, the ionising cross-section of the gas drops, leaving the 
large fraction of the gas at these distances neutral. The ionisation 
potential bins used in Fig.[5]are, however, crude, and especially 
in the case of gas with IP > 200 eV (i.e. iV u ition)> finer binning 
could produce a less pronounced dip in the ionisation profile at 
20 eV-50 eV (Ny) with an extended tail beyond 200 eV. In the 
context of the step-function-like density profile described above, 
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a more extended, shallower density profile would by produced 
by a smoother density profile. 

A more detailed investigation into the environmental condi- 
tions required to produce the observed ionisation profiles shown 
in Fig. [5] would require detailed modelling of the environment 
and radiation field produced by the GRB. Our ability to carry 
out such analysis at this time is, however, limited by the cur- 
rent lack of understanding of the conditions of the gas close to 
the GRB. Most of the lines measured in optical spectra trace 
material far from the GRB site, leaving us with little insight on 
the ionising flux of the radiation field, the structure and density 
of the environment, nor its temperature and metallicity. Greater 
constraint on the environmental properties within a few tens of 
parsecs to the GRB are, therefore, needed before an attempt to 
model the ionisation profile along the line-of-sight to the GRB 
can be made. 

5. Summary and Conclusions 

Using optical and X-ray spectra together with optical photomet- 
ric afterglow data, we have combined the results from both spec- 
troscopic and broadband SED analysis for a sample of 29 GRBs. 
This has allowed us to probe the properties of the host galaxy in- 
terstellar medium much deeper than would be possible with the 
analysis of optical or X-ray data alone. 

From a sample of 18 GRBs we confirm the decreasing trend 
in the iron-to-zinc relativ e abundance as a function of log iVz n ii 
previously observed (e.g. Savagliol l2006l) . and further find this 
correlation to be supported by independent dust extinction mea- 
surements from GRB afterglow SEDs. 

In comparing UV absorption lines of singly ionised metals, 
which are a tracer of neutral gas, to the soft X-ray absorption in 
the afterglows of 26 GRBs, which is representative of the total 
gas along the line-of-sight, we found the latter to be significantly 
larger. After taking careful consideration of systematic effects, 
we found the most likely interpretation for these observations to 
be that ~ 90% of the host galaxy gas along the line-of-sight to 
the GRB is ionised. 

Evidence for significant column densities of ionised gas 
wi thin GRB host galaxies has already been pre sented, such as 
in iProchaska et al) (l2008l) and iFox et alJ d2008l) . In both these 
papers, the velocity profiles of strong N V absorption in the host 
galaxies of a handful of GRBs were found to be consistent with 
arising from highly ionised circumburst gas photoionised by the 
GRB. However, from our analysis we found that only a negli- 
gible fraction (~ 0.01 %) of the X-ray absorbing gas is traced 
by th e highly ionised metals measured in Prochaska et al. (2008) 
and lFox et al.l d2008l) . This leads us to postulate that the majority 
of the ionised gas probed by our soft X-ray measurements is in 
an ultra-ionised stat e. 

IProchaska et ail |2008) estimated that a cold gas (T ss 
10 4 K), with constant density profile would be mostly ionised 
by the GRB to N V out to a radius of r w 10 pc. The indica- 
tion of our results that N V only traces a small fraction of the 
ionised gas within the GRB host galaxy thus suggests that most 
of the ionised gas lies inside the region where N V is detected, 
where the gas would experience a greater level of photoioni- 
satio n. The radius and d ensity of the ionised cloud estimated 
by IProchaska etaD {2008) cannot account for the large column 
densities of ionised gas suggested by our X-ray measurements. 
However, both sets of observations may be satisfied by imposing 
a wind-density profile ( n oc r 2 ), such that a larger column of gas 
could be photoionised to very high states within a few parsecs of 
the GRB, while still maintaining a significant column of N V. 



The recent surge in early time photometric and spectroscopic 
data available in the Swift era of GRB observations has pro- 
vided a significant better handle on the conditions of the host 
galaxy material along the line-of-sight to a GRB, which shall 
be well complimented by future observations with XShooter 
at VLT. Observations of variability in the fine-structure lines 
of singly ionised elements, such as Si II and Fe II are plac- 
ing limits on the distances of the n e utral absorbing gas (e.g 
IVreeswiik etafl 120071: iLedoux et all 120091) . and th e detailed 
analys is of strong-ion absorption lines presented in IFox et al] 
(20081) and IProchaska et all d2008l) provide some of the first at- 
tempts at investigating the GRB circumburst environment. 

From the broadband photometric and spectroscopic data pre- 
sented in this paper we have been able study the relation between 
the neutral and ionised gas within GRB host galaxies, providing 
greater insight into the affect a GRB has on its surrounding en- 
vironment, and on the importance of the numerous absorption 
components on our observations. The lack of optical absorption 
lines at lower IPs leads us to conclude that either the large frac- 
tion of the ionised gas along the line-of-sight to the GRB is in 
an ultra-ionised state, or that further consideration on our under- 
standing of the absorption of the GRB optical and X-ray after- 
glow is needed. An ultra-ionised gas would produce a signature 
in the soft X-ray energy bands, where absorption from oxygen 
and carbon is already prominent. It will, however, require future, 
high sensitivity X-ray spectroscopic missions such as Xenia to 
verify the presence and abundance of such highly ionised gas. 
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